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ABSTRACT 

It was found by Amati et al. in 2002 that for a small sample of 9 gamma-ray 
bursts, more distant events appear to be systematically harder in the soft gamma- 
ray band. Here, we have collected a larger sample of 65 gamma-ray bursts, whose 
time integrated spectra are well established and can be well fitted with the so 
called Band function. It is confirmed that a correlation between the redshifts (z) 
and the low-energy indices (a) of the Band function does exist, though it is a 
bit more scattered than the result of Amati et al. This correlation can not be 
simply attributed to the effect of photon reddening. Furthermore, correlations 
between a and E^eak (the peak energy in the uF^, spectrum in the rest frame), a 
and £^iso (the isotropic energy release), a and L^so (the isotropic luminosity) are 
also found, which indicate that these parameters are somehow connected. The 
results may give useful constraints on the physics of gamma-ray bursts. 

Subject headings: Gamma-ray burst: general — X-rays: bursts — Methods: statistical 
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Introduction 



Correlation analysis plays an important role i n leading us understanding the physics 
of astronomical processes (Dyson & Schaefer Il998l ). The correlations between the spectral 



parameters of the prompt emission of gamma-ray bursts (GRBs) and other related 
parameters help to reveal their nature. In this aspect, several empirical relationships 
have been found previously. The most famous one is the correlation between the rest 
frame peak energy (-Epeak, note that we will use E°^l^^ to designate the corresponding 
peak energy in the observer's frame hereinafter) of the uF^ spectrum and the bolometric 
iso tropic energy release (-Eiso) during the burst, also known as the Amati relation (Amati et 



al. 

col 

al. 



2002 



Amati 



20031 ). Significant correlation is also found to exist between -Epeak and the 



imat ion-corrected ener gy, wh ich is usually called the Ghirlanda relation (Ghirlanda et 



2004 ). Liang & Zhang ( l2005l ) further established a three-parameter correlation between 



Epeak, Eiso and the jet-break time in the afterglow hght curve. The fourth relationship 
that should be mentioned is between -Epeak and t he iso tropic peak luminosity (Liso), i.e.. 



2004 ). Very recently, for a sub-sample 



the so called Yonetoku relation (Yonetoku et al. 
of long GRBs with known redshifts and with a plateau phase in the afterglow, a tight 
three-parameter correlation has also been found among the end time of the plateau phase 



(in the GRB rest fra me), t 



release (Xu & Huang 120121 ). 



le corresponding X-ray luminosity and the isotropic 7-ray energy 



However, the nature of many of these relationships has been under debate since their 



discovery. Some be 
the burst (Schaefer 



ieve t hat they follow from the physics of the emission process durin g 

Ghirlanda et al. 



2004 : Bosnjak et al. 



2008; Amati et al. 



2009 



2010f ). 



Especially, a few of the correlations (such as the Yonetoku relation) are found to exist in 
time-resolved analysises of individual GRBs. The time-resolved behaviors are largely similar 
for different GRBs and are also consistent with the time-integrated relations (Ghirlanda 
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et al. I2OIOI : Ghirlanda et al. 1201 ll : Nava et al. |2012[ ). An intrinsic physical origin for these 
correlations is thus strongly hinted, since it is quite unlikely that instrumental selection 
effects or other observational biases can play the major role in a single burst. However, it 
should be noted that observational selec tion effects still can not becompletely excluded yet 



(Nakar & Piran 



2005 



; Butler et al. 



2007: Shahmoradi & Nemiroff 



2OI1I) 



The origin of the prompt em ission of GRBs is still an open question. Within the 



standard fireball framework (Piran |200j), many mechanisms (see Zhang & Meszaros 2002 for 
a review) have been proposed based on different ingredients of fireballs (e.g. kinetic energy 
dominated or Poynting flux dominated) and different emission processes (e.g. synchrotron 
emission from internal shocks, photospheric emission etc.). It has been suggested that both 
-E'peak and Eiso should be somehow linked to the bulk Lorentz factor (P) of the fireball in 
most scenarios (Amat i 200 6). This can potentially give an interpretation to the Amati 
relation. Liang et al. (120 lOf ) found a tight Pq - -^7,150 correlation by using the initial Lorentz 
factors constrained from afterglow analysises. They suggested that the correlation may be 
useful to pin down the physics of GRBs in the prompt phase. For example it may reflect 



the angular structures of the GRB jets. Recently Ghirlanda et al. (j2012al ) also connected 
the bulk Lorentz factors with observable parameters (-Epeak, Eiso, Lisa) and commendably 
explained the -Epeak - -Eiso, -f^peak - Liso correlations. 

The above correlations have prompted GRBs as potential "stan dard c andles" in the 



universe. Many studies have been devoted to this topic. Liang et al. ( 120081 ) used the SN la 



as a first-order standard candle to calibrate the GRB correlations, trying to avoid the so 

3e co smology depe ndent 



called circularity problem (the calibration of GRB correlations may 
due to the lack of a sufficient low-redshift GRB sample). Dai et al. ( 12004| ) and Amati (j2006l ) 
have used the Amati relation to constrain the cosmological models and obtained results 
consistent with those from the SNe la method. While the reliability of the Amati relation 
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for cosmological use is still under debate (Collazzi et al. |2012[ ). Schaefer (120071 ) argued 
that the consistency between GRB cosmology and conventional cosmology is robust, and 



that further improvements could be expected. Actually, Wright ( 120071 ) have used the GRB 
data as a supplement to the conventional data of the supernovae, acoustic oscillations, 
nucleosynthesis, large-scale structures, and the Hubble constant, to constrain the dark 
energy property. He concluded that the GRB Hubble diagram does help to break the 
degeneracy between the case of w 7^ 1 and the case of Qk 7^ 0. However, the large scatter 
of current GRB relations is still a serious problem for their applications in cosmology. 
Ti ghter or completel y new correlatio ns ar e still the goal ofmany researchers (Yonetoku et 



al. 



2010 



et al. 



Wang et al. 



2011 



20m . 



Qi & Lu 



2012 



Xu & Huang 



2012 



Mangano et al. 



2OI2I : Zhang 



We note from a plot in Amati et al.(2002) that an obvious correlation seems to exist 
between a (the low-energy spectra index of the time integrated GRB spectra fitted with 
the Band function) and the redshift z. We call this correlation the a - z relation. If the a - 
z relation really holds, then one can derive the redshifts of GRBs simply from the observed 



spectrums. It would be a new anc 
in many further studies (Guidorzi 



useful method to derive redshifts, which can be applied 



2005 



20081 ). However, the figure in Amati 



Curran et al. 

et al. only contains nine GRBs. No matter whether the a - z relation is intrinsic or not, it 
deserves to be tested by a larger sample. Today, the number of GRBs with known redshifts 
has been significantly increased thanks to extensive follow-up observations. In this study, 
we will examine the a - z relation by using a much expanded sample. We will also explore 
the correlations between a and other parameters, such as -Epeak? -^'iso? and Liso- 

Our paper is organized as follows. In Section 2, we describe the sample used for the 
study. In Section 3, the correlations among a variety of parameters are explored. It is found 
that there is indeed a correlation between a and z, although it is a bit more scattered 
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than the result of Amati et al. Correlations between a and -Epeak, -E'iso, -Z^iso are also found. 
Theoretical implications of these correlations are then explored in Section 4. Finally, 
Section 5 presents our main conclusions. Possible applications of the a - z relation are also 
discussed. 



2. Sample 



To fit the time-averaged spectra of GRBs, people usually use three kinds of functions: 
the si mple power-law function, the cutoff power-law function and the Band function (Band 



et al. 



19931 ). The simple power-law function is 



f{E) oc E- 



(1) 



where E is the photon energy and a' is the power-law index. The cutoff power-law function 
is 



/(E)(xE-"'exp(-^^^— ^). 



E. 



(2) 



peak 



The Band function is usually expressed as 



exp( 



-g(2-a) - 



E < 



-(a-/3)Epcak N 

2-a 



fiE) oc { 



(3) 



where a is the power-law index in the low-energy range, and /3 is the power-law index in 
the high-energy range. Note that generally a and f3 are positive in our notation. 

Each function has its advantage in specific spectrum fit. The derived spectral 
parameters based on different spectral functions are also systematically different. Spectrums 
of GRBs in the BATSE catalog are often well fitted by the Band function, while the spectra 
of GRBs observed by Swift are u sually best fitted by the power-law or the cutoff power-law 
functions (Sakamoto et al. 1201 if ). The difference may be caused by the relatively narrow 
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energy response of the Swift/BAT detector (15-150 keV). Band et al. (1993) have shown 
that even when the intrinsic spectrum of a burst is a Band function, the observed spectrum 
might still be well fitted with a simple power-law function or a cutoff power-law function 
if ^'peak is outside the energy range of the detector. This is intelligible because there is no 
suffici ent d ata on both sides of E°^l-^ to credibly constrain a function with a break (Krimm 



et al. 



20091). So, in order to obtain a good spectrum of a GRB, the detector should have 



a wide energy response. A wide energy band also help to get a good fit result with the 



Band function. Note that other factors, such as the sens itivity of the detector (Band 



20031) 



20091 ). may also affect the correct 



and the signal-to-noise ratio of the bursts (Krimm et al. 
derivation of spectral parameters. A less sensitive instrument or a low signal-to-noise ratio 
of the burst will lead to data with bad quality, which could not be used to credibly derive 
the spectral parameters. In this study, since we are mainly concentrating on the spectrum 
features of GRBs, we will only select those GRBs with high spectrum quality, i.e., they 
should be observed by a detector with wide energy response and high sensitivity so that the 
Band function parameters can be well constrained from the observations. Additionally, the 
redshifts of these GRBs also need to be measured. 

As discussed above, we need the parameters of the Band function derived from 
observations for the purpose of our study. In principle, the sample should be composed 
of GRBs with known spectroscopic redshifts and high-quality spectral data so that their 
spectrums can be well fitted with the Band function. Many GRBs observed by Swift 
have redshift data, but unfortunately the narrow energy band of the Swift/BAT detector 
(15-150 keV) seriously limits their usage in our study. As a result, our sample only contains 
a few GRBs from the Swift catalog. They were generally observed simultaneously by 
Swift/BAT and other wide-band detectors such as the Konus/Wind or Fermi/GBM, so that 
high-quality spectral data are available. 
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With the selection criteria mentioned above, we have collected 65 GBRs. The redshifts 
and spectrum parameters are listed in Table 1. For these GRBs, the spectrum are all 
better fitted with the Band function than with other functions according to the Ax^ 
cri terion (an improvemen t of 6 units in for a change of 1 degree of freedom, Band et 

20111). In fact, all the GRBs in our sample have been frequently 



al. 



1993: Sakamoto et al. 



us ed by other author s to test vari ous o ther correlatio ns (G hirlanda et al. 



2004 



Nava et 



al. 



2008 



2011b 



2012: Kann et al. 



2OI2I ). The spectrum parameters of 



20 Id : Tsutsui et al. 

the same GRB are generally consistent when the event appears in more than one article 
mentioned above. 

The distributions of ^'peak ^^^1 a of our sample are shown in Figure 1. The mean value 
of £^pgak is 282 keV and the mean value of a is 1.07 (with a standard deviation of o = 0.3 24). 



These values are similar to t 
Schaefer 



2003 



Kaneko et al. 



2000 



lose derived from the BATSE GRBs (Preece et al. 
20061 ). Note that the value of -Epcfk is typically well outside 



the energy band of the Swift/BAT detector, which clearly demonstrates the inaptness of 
most Swift GRBs for the purpose of our study. In Table 1, i^iso and Liso are calculated in 
the energy range of 1 - 10^ keV in the GRB rest frame (a scheme proposed by Amati et al. 



2002), and k-correction has been applied in the integration. GRBs are usual 
long (T90 > 2s) and short-duration (T90 < 2s) categories (Kouveliotou et al. 



y cla ssified as 



I993I). Under 



this classification scheme, all our events are "long GRBs" except for GRB 090510. 

For our sample, the observational data are from BeppoSAX, Konus/Wind, HETE-2 
and Fermi/GBM, which have broad energy coverage. The sample are heterogeneous in 
terms of the different instrumen ts, th us we need to consider the possible effects of the 



heterogeneity. Ghirlanda et al. ( l2008l ) analyzed the distributions of GRBs detected by 
different instruments (Swift, BATSE, HETE-2, Konus/Wind and BeppoSAX) in the 
-^peak ^ fluence plane. They found that the distribution of the heterogeneous sample is not 
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seriously affected by the trigger sensitivity (the minimum flux to trigger a burst). Except 
for the Swift sample (luckily not included in our sample), the distribution is also not 
affected by different spectral threshold (the minimum fluence required to get the spectrum 
and constrain the peak energy). For Fermi/GBM and BATSE long GRBs, they have similar 
distributions of fluence, -Epeak' peak flux. It is also found that the Fermi/GBM bursts 



generally have harder low-energy spectral indices 



but the difference is very slight (Nava et al. 



a) with respect to the BATSE GRBs, 



2011al ). So, although the heterogeneity of the 



bursts in our sample might do lead to some complex bias, it should not be too serious in 
general. 



3. Correlations 

We first investigate the a - z relation mentioned by Amati et al.(2002) with our larger 
sample. The original a - z relation in previous article is log a = (—0.78 ± 0.13) log(l + z) + 
(0.39 ± 0.04). In Figure 2, we plot our sample in the log a — log(l + z) plane. It confirms 
the decreasing trend of a with z, though the data points seem dispersive. Also, the slope 
of this correlation becomes flatter when the sample is expanded. The best fit result of the 
correlation is now 

log a = (-0.42 ± 0.07) log(l + z) + (0.11 ± 0.02). (4) 

The correlation coefficient is = 0.35. The associated p-value (the probability that such a 
correlation is simply formed by chance) is 2.3 x 10^''. In Figure 3, we plot our sample in the 
log (3 — log(l + z) and a — f3 planes. Massaro et al. once suggested that the a - z relation 
could be a manifestation of the selection effect: brighter GRBs, which can potentially be 
detected at higher redshifts, should be softened due to the e xpansion of the Universe, thus 



they might have a flatter photon index (Massaro et al. 



20021 ). If this explanation is correct. 



then both a and /3 should have similar correlation features with z, and additionally, a and 
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/3 themselves should be correlated. However, Figure 3 shows clearly that no correlation 
exists between /3 and z. Further more, the right panel of Figure 3 also shows that there is 
no positive correlation between a and /3. 

In fact, from theoretical aspect, the power-law spectral indices of the Band function 
should not be affected by cosmological redshifting. Assuming that the spectrums at rest 
frame and observer frame have similar function form as in Eq. (3), and denoting the 
spectral parameters in the rest frame and observe frame with a subscript "rest" and "obs" , 
then according to the conservation of photons in a given unit energy interval, we have 

diw^ . (5) 

dE dE 4 + ^ ^ 

This equation can be re-written as 



exp( « ( 1 )i-.o.^-... exp( (6) 

-^peak J- -|- ^ -^peak 



where -Epeak = (1 + ^)-^pcak- -^l- (6) leads to Orest = «obs- It indicates that the cosmological 
redshifting does not modify the value of the power-law index. So the a - z correlation 
should not be due to the reddening of photons. Actually, on the contrary, the a - z relation 
revealed in Figure 2 indicates that the low-energy spectrum gets harder with the increasing 
z. Additionally, the facts that /3 does not correlate similarly with and that /3 does not 
positively correlate with a (Figure 3) also strongly indicate that the a - z correlation 
should be due to some intrinsic mechanism. 

However, it should be noted that if -Epeak is shifted to be near the low-energy limit of 
the passband, then the low-energy data may not be sufficient enough to give the correct 
spectrum index. This kind of bias could be presented in the value of a. For example, a 
simple extrapolation of the Band spectra from the Fermi/GBM band to the Fermi/LAT 



band would systematically over-predict the observed flux (Ackermann et al. 



20121 ). This 



can be explained by the softer /3 values or intrinsic spectral breaks at energies > 40MeV, a 
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similar effect as discussed here. In Figure 4, we show our sample on the -Epefk ^ (1 + -2) plane 
to see to what extent this kind of bias might affect our results. Luckily, for most detectors, 
the lower limits are below 10 keV. As a result, actually all the GRBs lie significantly above 
the corresponding detector limit. So our results would not be seriously affected by this bias. 



Sakamoto et al. (120091 ) have found an empirical relation between the power-law 
spectum indices of the Swift/BAT GRBs (F^^^) and Epeak: log^peak « In Figure 5, 

a vs -Epeak and (3 vs -Epeak are plotted for the GRBs in our sample. We find that an obvious 
correlation exists between a and -Epeak, which can be best fitted as, 

log a = (-0.14 ±0.03) log Epeak(keV) + (0.33 ±0.07). (7) 

The corresponding correlation coefficient is = 0.30 and the p-value is 4.3 x 10^^. We 
call this relation the a - -Epeak relation. It implies that a GRB with a larger -Epeak tends 
to have a harder low-energy spectrum. However, we should notice that -Epeak itself is 
mathematically related to alpha through the definition of the e-folding parameter, Eq, by 
Epeak = -Eq * (2 — a). Therefore, the a - -Epeak relation may be partially due to the definition. 
No correlation exists between /3 and -Epeak? as shown in the right panel of Figure 5. 

We have studied the correlation between the spectral parameters and the total radiated 
energy. In Figure 6, we plot the sample in the log a - logE^iso plane. As clearly shown in 
the left panel of Figure 6, we see evidence of a correlation between i^iso and a (called the a 
- E[so relation hereafter). The best fit result is 

log a = (-0.08 ±0.01) log Ei,o(erg) ± (4.3 ±0.55). (8) 

The corresponding correlation coefficient is = 0.53 and the p-value is 3.8 x 10^^°. It 
implies that the more energetic one GRB be, the harder its low-energy spectrum is. But 
interestingly, for the high-energy index (/3), no trend of correlating with E^so can be found 
(see the right panel of Figure 6). 
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We have also investigated the correlation between the spectral parameters and the 
isotropic luminosity. The results are plotted in Figure 7. The left panel of Figure 7 shows 
that a correlates with Liso tightly. The best fitted function is 

loga = (-0.1 ± 0.01) logLiso(erg/s) + (5.2 ±0.6). (9) 

The corresponding correlation coefficient is = 0.63 and the p-value is 3.8 x 10"^^. Such 
an a - Liso relation is even tighter than the a - E\^o relation displayed in Figure 6. On the 
other hand, no correlation can be found between /3 and Liso, as shown in the right panel of 
Figure 7. 

Using our expanded sample, we can also examine the Amati relation. The results are 
plotted in Figure 8. From the left panel, we see that -Epcak and E-^^a are highly correlated. 
Our best fit result is: 

logEpeak(keV) = (0.46 ± 0.04) logEiso(erg) - (22 ± 2.2). (10) 

The corresponding correlation coefficient is = 0.71 and the p-value is 2.6 x 10^^^. The 
slope is ~ 0.5, consistent with previous studies. In the right panel of Figure 8, we plot the 
diagram of L\^o vs 2. It is shown that there is a tendency that the luminosity increases with 
increasing z. This could be due to the selection effect that faint GRBs are difficult to be 
detected at higher redshifts. The evolution of Ljso with respect to z can be fitted with a 
power-law function of 

log Liso(erg/s) = (2.56 ± 0.47) log(l + z) + (51.26 ± 0.14), (11) 

for which the correlation coefficient is = 0.39 and the p-value is 2.3 x 10^^. 

For our sample, there is also a clear correlation between -Epeak and z. Figure 9 shows 
that for more distant GRBs, the -Epeak value is generally higher. But generally speak ing, t he 
correlation is not so tight, roughly consistent with previous studies (Mallozzi et al. 



Wei & Gao 



20031) 



1995 
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As a comparison, we have also selected all the Swift GRsj^ with known redshifts, 
trying to examine any possible correlations between the spectral parameters and the 
redshifts. For the majority of these Swift GRBs, the spectral parameters of the Band 
function are not available due to the narrow energy range of the Swift/BAT detector. Their 
spectrums can only be fitted with a power-law or a cutoff power-law function. For these 
events, the derived power-law spectral indices may somewhat be similar to the lower-energy 
spectral indices of the Band function. However, somewhat unexpectedly, we find that there 
is no correlation between the power-law index and the redshift. It may indicate that the 
narrow energy range of the Swift/BAT detector has seriously limited our determination of 
the spectrum, even at the low-energy segment. 



4. Physical Implications 



For t he syn chrotron emission from shocked electrons in magnetic fields (Rybicki & 



Lightman 



19791 ). the low-energy spectrum index a should range from 2/3, in the case of 



optically thin synchrotron emission (Katz 19941 ). to 3/2, wh en th e distribution of electrons 



is characterized by a power-law index of —2 (Preece et al. Il998l ). However, a notable 
fraction of the observed GRB spectra are harder than the optically thin sync hrotron spectra 



(i.e. with a < 2/3). This is the so-called "death line" problem (Preece 



19981) 



Thermal emission may provide a solution to the "death line" prob lem. Ryde (120051 ) 
analyzed a sample of 25 pulses in the catalog of Kocevski et al. (120031 ). He found that 
thermal emission, combined with a nonthermal component, is ubiquitous for GRBs of 
various spectral shapes and timing characters. So, the GRB radiation may consist of 
both blackbody emission from photosphere and nonthermal component from relativistic 



^ http : / / swift . gsf c . nasa . gov / do cs/ swift / archive / grb _t able / 
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outflows (Ryde 2004; Fan et al. |2012| ). The measured value of a then should depend on 
the relative strength of the thermal component and the nonthermal component. Since the 
Rayleigh- Jeans portion of a blackbody spectrum is a power law function with a = —1, it 
may lead a to exceed the synchrotron range. For example, if the thermal component is 
strong enough to become the dominant component, then the resulted spectrum will have a 
small a. 



Strong photospheric emission at gamr na-ray wavelengt hs is 



scenarios. In u sual outflow models (Piran 



1999 



Meszaros 



predict ed in most GRB 



2002; Piran 



2004; Rees & 



Meszaros 



20051 ). the typical temperature of the photosphere is T oc L^/^. It indicates that a 



higher luminosity usually corresponds to a stronger thermal component, which then leads 
to a smaller a. This may naturally give an explanation to the a - Liso correlation, as 
displayed in Figure 7. At the same time, since both Liso and Eiso are indicative of the 
strength of GRBs, it is easy to understand that these two parameters should be positively 
correlated. Combining the a - L^so correlation and the Liso ^ -^'iso correlation, we can then 
easily explain the a - Eiso correlation. 

In Figure 10, we illustrate the relationship of all the correlations studied in our article. 
On the left side, from the a - E'iso correlation (Figure 6) and the Amati relation (i.e. the 
i^peak - Eiso correlation. Figure 8), we can get the a - Ep^ak correlation (see Figure 5, of 
course, the effect of deflnition as mentioned before cannot be ignored either). On the right 
side, from the a - L^so correlation (Figure 7) and the Liso ^ z correlation (a selection effect, 
as discussed in the previous section, also see Figure 8), we can get the a - z correlation 
(Figure 2). Finally, from the a - -Epeak and a - z correlations, we then have the L^peak ^ z 
correlation (Figure 9). 

However, we must bear in mind that the above analysis is only qualitative. When we 
check the correlations in detail, we get some difference. For example, Eq. (9) (i.e. the flt 
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result of Figure 7) indicates that logo; oc — O.llogLiso, and Eq. (11) (i.e. the fit result of the 
left panel of Figure 8) indicates logLiso oc 2.561og(l + z). Combining these two equations, 
we then get logo; oc — 0.261og(l + z). However, such a result is slightly different from the 
best fit result of Figure 2, which actually gives logo; oc — 0.421og(l + z). The difference 
indicates that there might be other factors that are playing their roles. They need to be 
studied further in the future. 

The largest uniform sample of GRBs with measured low-energy spectral indices (i.e. 
the a parameter of the Band function) should be the BATSE GRBs. However, the redshifts 
of most of the BATSE GRBs are unknown. Fortunately, we could estimate the redshifts 
from some statistics relations. We have calculated the pseudo redshifts of the BATSE 
GRBs by using the so called Yonetoku relation (the -Epeak - -^iso correlation, Yonetoku et 
al. 2004). We then can plot the BATSE GRBs on the logo; - log(l + z) plane. The results 
are shown in Figure 11. In this figure, we have also plotted our sample as the red dots. 
It is interesting to see that the BATSE GRBs also show a correlation between a and z, 
i.e. logo; oc — 0.231og(l + z). The behavior of BATSE GRBs and our sample are generally 
similar, though the slope is slightly different. It indicates that the a - z correlation studied 
here may also hopefully act as a new redshift indicator. 



5. Discussion and Conclusions 

We have collected a large sample of 65 GRBs from the literature, of which both the 
redshifts and the Band function spectral parameters are available. Using this greatly 
expanded sample, we investigated the correlation between a and z, as first proposed by 
Amati in 2002. It is confirmed that the a - z correlation does exist, although it seems 
more scattered in our case. At the same time, it is interesting to note that the high-energy 
spectral index, (3, does not correlate with z. A few other correlations are also found to exist 
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for our sample, including the ^^peak - -Eiso, ol - E-^, a - Li^o, Li^o - z, a - -Bpeak, -£^peak 
- z correlations, etc. Among all these correlations, the a - L^so relation is the tightest. 
Again, although the low-energy spectral index (a) shows an obvious correlation with many 
parameters such as £^iso, -^iso, z, £^peak, the high-energy spectral index (^) does not correlate 
with any of these parameters. 

The a - z correlation indicates that the "death hne" problem is more serious at 
high redshifts. We have shown that this correlation can not be simply attributed to the 
photon reddening induced by cosmological expansion. Amati (2002) suggested that the 
correlation may hint that the radiative cooling occurs more actively in GRBs at smaller 
redshift. However, the intrinsic factor that leads to such a difference of cooling is completely 
uncertain. According to the statistics of Ryde, thermal emission may be a common 
component in GRBs. Here we prefer to use the thermal component assumption to explain 
the correlation. For GRBs with a higher luminosity (i.e., a larger Ligo), the relative strength 
of the thermal component with respect to the non-thermal component is also higher, leading 
to a smaller a. At the same time, at higher redshifts, only those GRBs with a relatively 
higher luminosity are likely to be detected (i.e., the Liso - z correlation, as displayed in 
Figure 8, which is actually a reasonable selection effect). Considering these two factors, the 
a - z correlation can then be naturally explained. Based on this idea, the other correlations 
mentioned above can all be easily understood (see the schematic illustration in Figure 10). 

For the BATSE GRBs, we have calculated their pseudo redshifts by using the Yonetoku 
relation (Yonetoku et al. 2004). When plotted on the logo; - log(l -|- z) diagram, it is found 
that the BATSE GRBs show similar tendency of correlation as our sample. This indicates 
that the a - z correlation might serve as a new redshift indicator in the future. When the 
Band function of the spectrum of a GRB is available, the low-energy spectral index might 
be used to give useful information on the redshift. 
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In many studies, when the spectrums of GRBs are involved, it is usually assumed 
that the intrinsic spectrum shape is a standard one, or at least it does not evolve with 



redshift. For example, Ghirlanda et al. (l2012bf ) simulated a population of bursts to inspect 



the impact of selection biases on the -Epeak - -Z^iso correlation of GRBs. They assumed 
that the GRB spectra has fixed low and high energy spectral indices over the range of 
redshift investigated. But according to our current study, the observed spectral parameter 
a actually evolves with redshift. This may either be due to the fact that the intrinsic a 
correlates with other parameters such Eiso, or that the intrinsic a itself evolves 

with redshift. In any case, there does not exist a standard shape for the intrinsic spectrum. 

GRBs can be approximately divided into two categories, short GRBs with the duration 
less than ~ 2 seconds and long GRBs with the duration larger than ~ 2 seconds. It is 
generally believed that short GRBs may be produced by the mergi ng of compact binaries. 



while long GRBs come from the collapse of massive stars (Woosley 



MacFadyen & Woosley 



19931 : Paczynski 



1998 : 



19991 ). It should be noted that most of the GRBs in our sample are 



long GRBs. The correlations investigated here thus are only applicable to long GRBs. It 
is an interesting problem that whether these correlations also exist for short GRBs. The 
problem should be solved when more and more short GRBs with measured redshifts and 
well determined spectrums are available. 



We thank the anonymous referee for valuable comments and suggestions. This work 
was supported by the National Basic Research Program of China (973 Program, Grant 
No. 2009CB824800) and the National Natural Science Foundation of China (Grant Nos. 
11033002 and J1210039). 
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Fig. 1. — Left panel: distribution of ^^pg^k of our sample. Right panel: distribution of the 
low-energy spectral index {a) of our sample. 
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Fig. 2. — Correlation between the low-energy spectral index and the redshift for our GRB 
sample. The solid line is our best fit. 
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Fig. 3. — Left panel: dependence of the high-energy spectral indice on redshift for the 
GRBs in our sample (GRB 081007 is not plotted since its error bar is too large). Right 

panel: a vs /3 for the same GRBs. 
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Fig. 4. — Distribution of the GRBs of our sample on the -Epeak ^ (1 + -2^) plane. The dotted 
hues correspond to the lower enegy limits of BATSE (yellow) and Konus/Wind (green) 
detectors. The limits of other detectors are not shown here because they are below 10 keV. 
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Fig. 5. — Left panel: ot vs -Epeak for the GRBs in our sample. The solid line corresponds to 
the best fit. Right panel: the corresponding /3 vs -Epeak diagram. 




Fig. 6. — Left panel: a vs E\^o for the GRBs in our sample. The solid line corresponds to 
the best fit. Right panel: the corresponding /3 vs i^igo diagram. 
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Fig. 7. — Left panel: a vs Liso for the GRBs in our sample. The solid line corresponds to 
the best fit. Right panel: the corresponding (3 vs Ljso diagram. 
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Fig. 8. — Left panel: the Amati relation (-Epeak vs -Eiso) for the GRBs in our sample. The 
solid line corresponds to the best fit. Right panel: Ljso vs z for the same sample. The solid 
line is our best fit. 
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Fig. 9. — Epeak (rest frame) vs z for the GRBs in our sample. The sohd hne corresponds to 
the best fit result. 
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Fig. 10. — A schematic illustration of the relationship of the correlations studied in our work. 
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Fig. 11. — A comparison of our sample (red dots) and GRBs from the BATSE catalog (black 
dots) on the logo; — log(l + z) plane. The pseudo redshifts of BATSE GRBs are derived 



from the Yonetoku relation (Yonetoku et al. 
our sample (red) and BATSE GRBs (black). 



20041 ). The solid lines are the fitting results of 
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Table 1. List of GRBs with known redshifts and Band function spectral parameters that 

define our sample. 



Name 


2 


a 


/3 


Range 
(keV) 


-'-^pcak 

(keV) 


(lOS^erg) 


T ■ ^ 

-^iso 

(lO^ierg/s) 


Mission'^ 


Ref. 


970228 


0.695 


1.54 ± 0.08 


2.5 ± 0.4 


2-700 


195 ± 64 


1.86 ± 0.14 


9.1 ± 2.18 


SAX/WFC 


1,2 


970508 


0.835 


1.71 ± 0.1 


2.2 ± 0.25 


2-700 


145 ± 43 


0.71 ± 0.15 


9.4 ± 1.25 


SAX/WFC 


1,2 


970828 


0.958 


0.7 ± 0.08 


2.1 ± 0.4 


30-10000 


583 ± 117 




25.1 ± 7.7 


RXTE/ASM 


2 


971214 


3.42 


0.76 ± 0.1 


2.7 ± 1.1 


2-700 


685 ± 133 


24.5 ± 2.8 


72.1 ± 13.3 


SAX/WFC 


1,2 


980326 


1 


1.23 ± 0.21 


2.48 ± 0.31 


2-700 


71 ± ,36 


0.56 ± 0.11 


3.47 ± 1 


SAX/WFC 


1,2 


980613 


1.096 


1.43 ± 0.24 


2.7 ± 0.6 


2-700 


194 ± 89 


0.68 ± 0.11 


2 it 0.67 


SAX/WFC 


1,2 


980703 


0.966 


1.31 ± 0.14 


2.39 ± 0.26 


50-300 


499 ± 100 




20.9 ± 4.86 


RXTE/ASM 


2 


990123 


1.6 


0.89 ± 0.08 


2.45 ± 0.97 


2-700 


2030 ± 161 


278.3 ± 31.5 


353 ± 123 


SAX/WFC 


1,2 


990506 


1.307 


1.37 ± 0.15 


2.15 ± 0.38 


50-300 


653 ± 130 




41.8 ± 13.3 


BATSE 


2 


990510 


1.619 


1.23 ± 0.05 


2.7 ± 0.4 


2-700 


423 ± 42 


20.6 ± 2.9 


61.2 ± 10.7 


SAX/WFC 


1,2 


990705 


0.843 


1.05 ± 0.21 


2.2 ± 0.1 


2-700 


348 ± 28 


21.2 ± 2.7 


16.5 ± 2.77 


SAX/WFC 


1,2 


990712 


0.433 


1.88 ± 0.07 


2.48 ± 0.56 


2-700 


93 ± 15 


0.78 ± 0.15 


0.746 ± 0.191 


SAX/WFC 


1,2 


991208 


0.706 


1.1 ± 0.4 


2.2 ± 0.4 


20-2000 


312.3 ± 5.1 


26.9 ± 3.7 


43.2 ± 3.8 


K/W 


2,3 


991216 


1.02 


1.23 ± 0.25 


2.18 ± 0.39 


20-2000 


763.6 ± 20.2 


85.1 ± 9.8 


113 ± 37.5 


BATSE 


2,3 


000131 


4.5 


0.69 ± 0.08 


2.07 ± 0.37 


50-300 


714 ± 142 


190.5 ± 30.7 


141 ± 55.9 


BATSE 


2,3 


000214 


0.42 


1.62 ± 0.13 


2.1 ± 


40-700 




0.8 ± 0.026 




SAX/WFC 


4 


000301C 


2.04 


1 ± 0.3 


2.3 ± 0.3 


20-2000 


1213 ± 303 


9.8 ± 0.9 




SAX/WFC 


3 


000911 


1.06 


1.11 ± 0.12 


2.32 ± 0.41 


15-8000 


1856 ± 371 


72 ± 15 


165 ± 28.9 


IPN 


2 


000926 


2.04 


1.1 ± 0.1 


2.43 ± 0.4 


20-2000 


306.9 ± 18.2 


33.1 ± 6.1 


47.3 ± 13 


K/W 


2,3 


010222 


1.477 


1.35 ± 0.19 


1.64 ± 0.02 


20-2000 


765.4 ± 29.7 


114.8 ± 7.9 


7.87 ± 0.45 


SAX/WFC 


2,3 


010921 


0.45 


1.49 ± 0.16 


2.3 ± 


30-700 


153 ± 31 


0.9 ± 0.1 


0.733 ± 0.133 


HETE-2 


2,4 


011121 


0.362 


1.1 ± 0.3 


2.1 ± 0.2 


20-2000 


1201 ± 126.7 


10.23 ± 1.6 




SAX/WFC 


3 


011211 


2.14 


0.84 ± 0.09 


2.3 ± 0.46 


40-700 


185 ± 22 


10.23 ± 1.4 


3.17 ± 0.,32 


SAX/WFC 


2,3 


020124 


3.198 


0.87 ± 0.17 


2.6 ± 0.65 


2-400 


390 ±113 


30.2 ± 3.6 


51.2 ± 20.3 


HETE-2 


2,4 


020813 


1.25 


0.94 ± 0.03 


1.57 ± 0.04 


20-2000 


710 ± 33.8 


87 ± 8 


25.8 ± 2.4 


HETE-2 


2,3 


020819B 


0.41 


0.9 ± 0.15 


2 ± 0.35 


25-100 


70 ± 21 




1.49 ± 0.323 


HETE-2 


2 


021211 


1.01 


0.85 ± 0.09 


2.37 ± 0.42 


2-400 


94 ± 19 


1.1 ± 0.13 


7.13 ± 0.99 


HETE-2 


2,4 


030226 


1.987 


0.89 ± 0.17 


2.3 ± 0.46 


20-2000 


349.5 ± 41.8 


14.8 ± 7.5 


8.52 ± 2.23 


HETE-2 


2,3 


030328 


1.52 


1.14 ± 0.03 


2.1 ± 0.3 


2-400 


328 ± 35 


36 ± 11.7 


11 ± 1.55 


HETE-2 


2,3 


030329 


0.169 


1.32 ± 0.02 


2.44 ± 0.08 


2-400 


79 ± 3 


1.74 ± 0.08 


1.91 ± 0.237 


HETE-2 


2,3 


030429 


2.66 


1.12 ± 0.25 


2.3 ± 0.46 


2-400 


128 ± 35 


1.78 ± 0.49 


7.6 ± 1.47 


HETE-2 


2,3 



-28- 



Table 1 — Continued 



Name 


z 


a 


/3 


Range 
(keV) 


.S'peak'' 
(keV) 


(lO^^erp) 


J ■ 

(10^1 erp/s) 


Mission'^ 


Ref. 


040924 


0.858 


1.17 ± 0.23 


2.3 ± 0.46 


20-500 


125 ± 11 


0.89 ± 0.2 


6.1 ± 1.1 


HETE-2 


2,3 


041006 


0.716 


1.37 ± 0.27 


2.3 ± 0.46 


30-400 


109 ± 22 


2.04 ± 0.6 


8.65 ± 1.36 


HETE-2 


2,3 


050401 


2.9 


0.99 ± 0.19 


2.51 ± 0.23 


20-2000 


458.2 ± 70.2 


42.7 ± 2.9 


203 ± 10 


K/W 


2,3 


050525A 


0.606 


1.02 ± 0.11 


3.12 ± 0.4 


15-10000 


130.1 ± 4.8 


2.4 ± 0.2 


7.23 ± 0.18 


K/W 


3,5 


050603 


2.821 


0.79 ± 0.06 


2.15 ± 0.09 


20-3000 


1334 ± 107 


61.7 ± 1.4 


2130 ± 220 


K/W 


2,3 


061007 


1.261 


0.75 ± 0.02 


2.79 ± 0.09 


20-10000 


965 ± 27 


101 ± 1.4 


109 ± 9.1 


K/W 


5 


061222A 


2.09 


1 ± 0.05 


2.32 ± 0.38 


20-10000 


1091 ± 167 


22.5 ± 0.94 


140 ± 38 


K/W 


5 


070125 


1.547 


1.1 ± 0.1 


2.08 ± 0.13 


20-10000 


934 ± 148 


95.5 ± 11 


324 ± 50 


K/W 


2,3 


070328 


0.372 


1.11 ± 0.04 


2.33 ± 0.24 


20-10000 


1052 ± 152 


1.95 ± 0.07 


2.77 ± 0.56 


K/W 


5 


071010B 


0.947 


1.25 ± 0.6 


2.65 ± 0.35 


20-1000 


101 ± 23 




6.4 ± 0.053 


K/W 


2 


080319B 


0.937 


0.86 ± 0.01 


3.59 ± 0.45 


20-7000 


1307 ± 443 


142 ± 3 


102 ± 9.4 


K/W 


5 


080721 


2.591 


0.96 ± 0.07 


2.42 ± 0.29 


20-7000 


1785 ± 223 


121 ± 10 


1038 ± 172 


K/W 


5 


080916C 


4.35 


0.91 ± 0.02 


2.08 ± 0.06 


10-1000000 


2268 ± 128 


309 ± 64 




Fermi 


3 


081007 


0.53 


2.1 ± 0.1 


10 ± 10 


25-900 








Fermi 


6 


081121 


2.512 


0.46 ± 0.08 


2.19 ± 0.07 


8-35000 


608 ± 42 


30.5 ±2.6 


195 ± 31 


Fermi 


5 


081222 


2.77 


0.9 ± 0.03 


2.33 ± 0.1 


8-35000 


630 ± 31 


25.2 ± 2.3 


94.9 ±3.1 


Fermi 


5 


090323 


3.57 


0.96 ± 0.12 


2.09 ± 0.19 


20-10000 


1901 ± 347 


398 ± 73 




K/W 


3 


090328 


0.735 


0.93 ± 0.02 


2.2 ± 0.1 


8-1000 


1133 ± 78 


19.5 ±0.9 




Fermi 


3 


090424 


0.544 


1.02 ± 0.01 


3.26 ± 0.18 


8-35000 


250 ± 3.4 


3.97 ±0.08 


11.2 ± 0.17 


Fermi 


5 


090510 


0.903 


0.8 ± 0.03 


2.6 ± 0.3 


8-40000 


4400 ± 400 






Fermi 


6 


090618 


0.54 


1.26 ± 0.04 


2.5 ± 0.25 


8-1000 


155.5 ± 11 






Fermi 


6 


090902B 


1.822 


0.61 ± 0.01 


3.87 ± 0.16 


8-1000 


798 ± 7 






Fermi 


6 


090926A 


2.1062 


0.693 ± 0.01 


2.34 ± 0.011 


8-1000 


268 ± 4 






Fermi 


6 


091003 


0.8969 


1.04 ±0.05 


3.58 ±1.21 


8-35000 


755.7 ± 70 


9.16 ± 2 




Fermi 


7 


091020 


1.71 


1.2 ± 0.06 


2.29 ± 0.18 


8-35000 


507 ± 68 


7.96 ±1.16 


32.7 ± 4.6 


Fermi 


5 


091127 


0.49 


1.25 ± 0.05 


2.22 ± 0.01 


8-35000 


51 ± 1.5 


1.61 ±0.03 


9.08 ± 0.22 


Fermi 


5 


100414A 


1.368 


0.37 ± 0.03 


3.74 ± 0.95 


8-35000 


1356.4±38.4 


76.2 ± 10 




Fermi 


7 


100621A 


0.542 


1.7 ± 0.13 


2.45 ± 0.15 


20-2000 


146 ± 23 


4.35 ±0.48 


3.17 ± 0.24 


K/W 


5 


100814A 


1.44 


0.74 ± 0.13 


2.73 ± 0.69 


8-35000 








Fermi 


7 


110213A 


1.46 


1.24 ± 0.17 


2.08 ± 0.05 


8-35000 


122 ± 20 


8 ± 0.4 




Fermi 


7 


110422 A 


1.77 


0.65 ± 0.036 


2.96 ± 0.12 


20-5000 


421 ± 8 


72 ± 0.3 




K/W 


7 
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Table 1 — Continued 



Name z q /3 Range ^'peak 

-Biso'" -f'iso'' Mission'^ Ref. 

(kcV) (kcV) (lO^^erg) (lO^^erg/s) 



110503A 1.613 0.98 ± 0.08 2.7 ± 0.3 20-5000 572 ±50 18 ± 1.4 181 ± 18 K/W 5 
110715A 0.82 1.23 ± 0.06 2.7 ± 0.3 20-5000 218 ±13 4.8 ± 0.2 K/W 7 

110731A 2.83 0.67 ± 0.15 2.64 ± 0.42 8-35000 1080 ± 133 45.5 ± 4 Fermi 7 



''Bpoak is the peak energy of uF,j spectrum in the GRB rest frame. 
°The integrating range of Eiso is 1 - 10^ keV in the GRB rest frame. 
'^The integrating range of Ljso is 1 - 10* keV in the GRB rest frame. 

'^The name of the mission from which the spectral parameters have been derived: SAX/WFC = BeppoSAX, K/W 
Konus/Wind, IPN = Interplanetary Network. 



Ref erence s. — (1) Amati et al. 1I2OO2I) : (2) Nava et a l ||200S|) : (3) Kann et al. ||2010|); (4) Ghirlanda et al. (^2004i 1: (5) Nava 
et al. 120121) : (6) Nava et al. J2011bl '): (7) Tsutsui et al. J2012r) . 
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